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1. INTRODUCTION 

Currently, a continuous growing demand for clean sustainable energy has pushed academic as well 
as industrial institutions to gain a deep understanding of solar cells. Much of today’s research focuses on the 
cells development, material selection and manufacturing with the aim to enhance performance and reduce 
cost. The current conventional photovoltaic is crystalline silicon known for its non-toxicity, abundance, high 
and stable cell efficiencies. It has been dominating the solar panels market due to its nature material 
technology and widespread knowledge, skills and infrastructure for manufacturing and assembly lines. The 
upper limit of crystalline silicon solar cells efficiency is reaching 26.7% [1]. However, they have lower 
energy conversion efficiency which make them less space efficient, limited flexibility and cannot be 
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manufactured in large scale. These limitations of silicon solar cells are some of the reasons that has led many 
researchers to explore other alternative materials to overcome these issues. 

Photovoltaic cells formed by organic polymers or organic molecules are solar cells potentially able 
to increase the energy efficiency compared to the first and second generations. They are well known for their 
lightweight, flexibility and low cost to fabricate where the manufacturing techniques vary from coating, 
spraying and printing [2]. On the contrary, the main disadvantages associated with PSC when compared to 
inorganic solar cells have low stability in the face of climatic variations and also a low energy efficiency. 
However, these factors had made the organic solar cells a challenging and popular field of research. The 
PSCs were able to achieve the highest efficiency of 2.5% as reported in [3], and continues to improve to a 
record level of 17.35% as presented in solar cell efficiency tables [1], [4]. One of the major obstacles that 
limit the energy conversion efficiency of organic solar cells lies in the intrinsic properties of the currently 
used active materials such as, conjugated polymers, preparation and control of their structure (distribution of 
different lengths of polymer chains in the material). In addition, the PSC tends to suffer from degradation 
when exposed to harsh ambient conditions (temperature, humidity and rain) as the solar cells are sensitive to 
high temperature [5], [6]. Several research groups are actively working to address these issues. For instance, 
a study to investigate the effect of intrinsic parameters which are influenced on the organic cells 
performance, such as Voc (open circuit voltage), Isc (short circuit current) and n (power conversion efficiency) 
was reported in [7]. Numerous models were also developed to represent the solar cell covering mainly the 
number of diodes, the shunt resistor and ideality factor. For example, the single diode model M1DR in [8] 
and two diodes model M2DR in [9]. We consider in the single diode model the following five parameters: 
the photocurrent Ipn, the reverse saturation current Io, the series resistance Rs, the shunt resistance Rsn, and the 
diode ideality factor n [10], [11]. The second problem is focused on the temperature effect of the (I-V) 
characteristic solar cell. This dependence is described by the sensitivity of the different electrical parameters 
of the organic photovoltaic generator (Voc, Isc, Nn, and Pmax) when the temperature varies. Indeed, the variation 
of the short circuit current in a remarkable manner, however the open circuit voltage is insensitive to the 
variation of the temperature [12]-[16]. In literature, several numerical methods were developed for extracting 
electrical parameters solar cells, such as the least square method, the Newton’s method, and the Levenberg 
Marquardt’s method [17]-[20]. As is justified in the literature, an accurate optimization method for 
parameters extracting is based on the minimization of root mean square error (RMSE) between the theoretical 
and experimental values to verify the effectiveness of the proposed model [21], [22]. 

The aim of the current study is to investigate the effect of temperature on the electrical parameters 
of the polymer solar cells (PSCs) based on poly (3-hexylthiophene) (P3HT) and [6,6]-phenyl C61-butyric 
acid methyl ester (PCBM). To further enhance the energy conversion efficiency, a developed model of the 
PSC device using single diode will be presented together with the proposed numerical optimization based on 
the combination of both the Least Square Method (LMS) and the Newton-Raphson (N-R) Algorithm. 


2. STRUCTURE AND OPERATION OF POLYMER SOLAR CELL 

The general structure of PSC device consists of a transparent substrate such as (glass, quartz, 
plastic), coated with an electrode, indium tin oxide (ITO), which should be transparent in order to allow the 
incoming photons to reach the photoactive layer. In addition, a bi-layer of two organic components, with one 
as the electron donor and the other as the electron acceptor, and finally one top metal electrode, as seen in 
Figure 1. This bilayer of PSC is the photoactive layer which composes of donor layer and acceptor layer. 

A schematic of the basic operation of a PSC is illustrated in Figure 2. the four main steps to produce 
electrical energy from solar energy by the solar cell are summarized as: 1) creation of electron-hole pairs so- 
called excitons by photons; 11) the excitons formed diffuse towards the donor-acceptor interface (D/A); 111) 
the minority charge carriers which are the holes in the donor region and the electrons in the acceptor region, 
move to the contacts for collection due to the internal electric potential/field in the junction; and iv) 
collecting The charge carriers associate at the electrodes and thus create an electric current. 

Figure 3 shows the operation mechanism in P3HT/PCBM cell. The exciton formed is generated, 
when the incident photon has energy greater than the bandgap energy E, of the P3HT. Thus, an electron in 
the highest occupied molecular orbit (HOMO) would be excited to the lowest unoccupied molecules orbital 
(LUMO), leaving a hole HOMO band. On the other hand, an efficient solar cell should have a wide 
absorption spectrum, to create as many pairs as possible of excitons charges due to the Shockley and 
Queisser limitation [22]. 
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Figure 3. A general operation mechanism in a bilayer P3HT: PCBM cell 


The created excitons can be separated if they encounter an internal electric field in the junction. 
Electrons move from LUMO donor to LUMO acceptor while the holes located in the HOMO acceptor and 
flows to HOMO donor. However, the diffusion length of charge carriers which is very short makes this 
process inefficient. Furthermore, the generated charge would be collected at electrodes and more efficient 
charge collection by utilizing metal oxide-organic interface to decrease surface resistivity and metal oxide- 
organic to facilite hole-extraction [23]-[25]. This transfer of carrier charges within the heterojunction, formed 
by the donor and acceptor of organic semiconductor D/A can produce an open circuit voltage V,.which is 
defined as reported in [26]: 


Woc = Estomo a E imo = Eoy (1) 


Where, EPomo and Ef;mo are the acceptor energy LUMO and the donor energy HOMO 
respectively, Eey 1s the exciton binding energy, and q is the elementary electric charge. 


3. MODELING OF PSC 
The one diode model of polymer solar cell under illumination and the equivalent circuit model is 
illustrated in Figure 4. The equation of current I as a function of voltage V is written as [20]: 





(2) 


q(V+IRs) 
I(V) = In- Io e ner - 1 -2 


Rsh 
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Figure 4. Equivalent circuit model of a polymer solar cell 


Where, I, and Io are the photocurrent and reverse saturation current, the series and parallel resistances are 
noted as R, and Rsp, the ideal factor symbol is noted n , the Boltzmann’s constant is noted k, and the 
temperature cell is noted T and expressed by (K). 

The photocurrent generated by the cell varies linearly on solar irradiation and is also influenced by 
the temperature according to the (3) [27]-[29]: 


lon = (Innes + Àsc(T = Le) 5) (3) 


Where, Ipnref is the photocurrent of the module at Standard Test Conditions (STC) (usually 25 °C and 
1000W/m7?), Ter is the reference temperature of the PSC which is given on the data sheet, Asc is the 
temperature coefficient of short current which is found on the data sheet and is expressed by (A/K), G is the 
irradiation on the device surface (in W/m7*), and G,e f is the reference irradiation. The relationship between 
lọ as a function of the temperature T is given by the empirical by Green [30]: 


I, = Aexp — (= ) (4) 


kT 


Where A is a constant equal to 1,5 x 108 mA. cm~?and E, represents the band gap energy. The expression of 


external quantum efficiency (EQE) is the ration between the electrical energy and the incident photon flux. 
Then, the EQE is given by [31], [32]: 





oN — Isc(A) 
EQE(%) = 1240 x Ox 100 (5) 


Where, ~,(T) is the monochromatic light intensity (mW /cm?). J,-(A) is the photocurrent density 
(mA/cm7), A is monochromatic wavelength of light (1m).The EQE equal to unit if all photons are observed 
by the semiconductor. The expression of incident photon flux which depends especialy with the body’s 
temperature T(K) is given by Planck’s law [33], [34]: 


2mhc2 1 


palT) = — a grkg (6) 


Where, h is the Planck's constant, kpg is the Boltzmann's constant, and c is the speed of light in vacuum, T is 
the temperature of the black body. The used Planck’s law is implemented in a MATLAB script which yields 
the photon flux of the black body, and it is represented in Figure 5. 

On the other hand, the generated current density J,. ,by the PSC with incoming solar radiation, can 
be calculated from the external quantum efficiency EQE(A,T). This evaluation is obtained by integrating 
over the entire solar radiation spectrum [35]: 


Ise = q Sy" EQE(A)p,(T) da (7) 


Where, ~,(T) is the photon flux of solar radiation spectrum under AM, 5. Then, the short circuit current can 
be deduced using the area of the illuminated dvice noted A. It s written as: Ise = A X J¢¢. 
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Figure 5. The incident photon flux for an AM) ,s5 solar radiation spectrum at 6000 K 


4. ELECTRICAL PARAMETERS OF AN PSC 

The most important parameter of a PSC is the power conversion efficiency(7), which represents the 
ability of the cell to convert light into electricity. It is defined as the ratio between the maximum power 
output Panay and the incident power Pine = 100 W/m? at (STC). It is written as: 


P FFXVocXI 
n — — x 100% — a SC 


unc unc 


x 100 (8) 


Where FF is the Fill Factor corresponding to the maximum power ratio on the ideal theoretical power. FF 
describes the part of energy loss in the form of Joule heat, which mainly determined by surface resistivity and 
it expressed as: 


FF = "a x 100% = “past x 100% (9) 


Voclsc oclsc 


The short circuit current J; is the current when the voltage across the device is zero, and the open-circuit 
voltage V,, is the voltage when the current in the device is zero. According to the (9) and (10), it is clear that 
a high efficiency could be reached with higher Ve, Ise and FF. The voltage V,. can be expressed as a function 
of the currents J, and Io. It can be written as [30]: 


Ve. = in (= = 1) (10) 


5. DEVELOPED METHOD OF EXTRACTING ELECTRICAL PARAMETERS 

Determining the final optimal parameters values requires first an estimation of the initial guess 
assigned to the extracted parameters 0) [ee Ip, Rs, Ryn; which affect the convergence of the Newton- 
Raphson’s algorithm (N-R). The initial parameters 0) are obtained from the current-voltage experimental 
curve using the least mean square (LMS) method. This technique is one of the most widely used methods for 
modeling experimental measurements by a predetermined analytical function given by (2). This method 
consists to minimize the mean absolute or relative difference between the N measures /m (Vi) and the 
calculated values of electric current J, (Vi) obtained by (2). The mathematical model of single diode for 
estimation five parameters in the form as: 


aE TEs) 1 V+IRs 


f Uli Va 0) = lpn — lo |e" mer " (11) 


The second step, we use the N-R algorithm to determine the final five unknown electrical 
parameters 0 |I, v lo Rs Rsh, n| which characterize the single diode model of the polymer solar cell. The N-R 
method is widely used in finding the root of nonlinear equations. This method uses the derivative of f(x) at x 
to estimate a new value of the root. The desired precision is reached by iteration. if 0;,, 1s an approximation 
to the root of f(0) = 0. The N-R is an iterative method that estimate the root of a real valued nonlinear 


function f(@) with an initial guess 05. This linearization process is obtained through the first order Taylor’s 
series expansion (linear in 0) to the first around the guess (6;), which is given below [36]: 
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f Gi41) = f(O) + (Oiza — ODF) + = 0 (12) 


Therefore, the next approximation to the root 6; in the ith iterative step is given by the iterative formula [37]: 


Oi41 = 0i — f (8:)/f' (Gi) (13) 


The procedure is repeated until |0;,, — 0;| < €, where £ is a perset tolerance. The N-R method is easily 
extended to the solution of set of n-nonlinear equations, where n > 2. Therefore, the parameters 0 are 
evaluated at each iteration according to the formula [21]: 


0i = bi- — DOi, x flo (14) 


Where J is the Jacobian matrix for the function f(@) given by the first-order contains the partial derivatives at 


each parameter of the vector 0;. It is defined by: Ji; = i A0 = 0; — 8;_, denotes the change in the value of 
j 


the vector of roots (step size). The used function @ to optimize the electrical parameters is the objective 
function ø, which is defined as the difference between both of the experimental /°*? and calculated [°* 
electric current values. The expression of @ can be given by: 


= ex he)? _ 
pll Vp 0) = Diy? — H") = XL, ef (0) (15) 
Where, e; = I; *P _ eal is the error between the experimental and theoretical electric currents ike “P and 6% 
respectively, and N is the number of measurements. To validate our adopted method and evaluates its degree 
of precision, we use the root mean square error (RMSE) which is written as: 


2 
RMSE = 2 N (e(l, Vp 9;)) (16) 


Where (l;, Vi) represents the ith measured value of N experimental data. In a MATLAB script, the flowchart 
of the developed method is shown in Figure 6. 
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Figure 6. Flowchart for estimation PSCs parameters using the developed method 
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6. RESULTS AND DISCUSSION 

In this section, we will present the results of the theoretical research related to the dependence of the 
PSCs energy conversion efficiency with the temperature and the intrinsic electrical parameters by using the 
single diode model. The simulated results of /-V characteristic, under the Standard Test Conditions (STC) at 
AM, s, were carried out using the developed method which is composed of two means. The first one is the 
least mean square method (LMS) so as to determine the initial electrical parameters 6) [Ion lo Rs, Rsn» n], 
and the second one is the Newton- Raphson’s algorithm (N-R) to evaluate the final optimal values. 

Figure 7 shows both the calculated and experimental curves of the PSHT: PCBM polymer solar cell. 
That is, the experimental curve was represented by the squares in the Figure 7, whose the data points (I;i, Vi) 
with (i = 1,...,36) are from the reference cell [38], what’s more, the calculated curve was schematized by 
the solid line. 

As it is shown in Figure 7, there is a very good agreement between the simulation results and the 
experimental data (li, Vi). In the developed method, the LMS method is used to determine the initial 
parameters so as to ensure the convergence of the Newton’s algorithm. The final values of the five 
parameters are substituted in (2), and are obtained while the error criterion A0 = |0,;,, — 0;| is minimized. 
Then the PSC’s parameters which are found by the developed method are listed in Table 1 and compared to 
the references work. 


@#eeeeeoeos e + Experimental 
t+ LSM Method | | 








0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Voltage (V) 


Figure 7. Calculated and experimental I-V characteristics of PSHT: PCBM polymer solar cell 


Table 1. Comparaison of estimated parameters for single diode model of PSC 
Electrical Parameters | Proposed method Ref[39] Ref[40] Ref[41] 


In (mA) 6.59 6.4 6.4 6.4 
1,(A) 1.46x10°!4 2.00x10°8  2.63x10 1° 2.21x10°" 
n 0.84 3.027 1.363 1.068 
Ryn(Q) 1038.13 36200 36300 36400 
R,(Q 38.85 107.42 35.46 38.24 


From the results presented in Table 1, it is clearly seen that the five electrical parameters evaluated 
by the proposed method are quite closer than those found in the references [39]-[41]. The highest precision of 
the extracted values is directly attached to the lowest value of the root mean squared error (RMSE), which 
measures the average of squared differences between prediction and real in a set of electric current values. 
The found values of the statistical error (RMSE) are presented in Table 2 and compared to the reference 
values [39]-[41]. 


Table 2. Statistical error (RMSE)of electrical parameters model 
Statistical Error | Proposed Method Ref[39]| Ref[40]| Ref[41] 
RMSE 1.0012x107 2.3423x10° 1.4936x10° 6.3944x107 


According to the statistical error, the obtained value by the proposed method is lower than those 
found by other works [39]-[41]. As a result, our proposed method is more accurate than the other methods. In 
addition, the program’s execution is fairly fast than the other programs. In what follows, we will study the 


Int J Pow Elec & Dri Syst, Vol. 12, No. 3, September 2021 : 1772 — 1783 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1779 


effect of the influence of many factors on the PSC operation, such as: temperature, series and shunt 
resistances, ideality factor and illumination. This effect will be shown through both the variation of I-V and 
P-V characteristics in order to optimize the energy conversion efficiency of polymer solar cell. 


6.1. Temperature cell effect on PSC parameters 

The Figure 8 shows the impact of the temperature cell ranging from 10°C to 90°C on the J-V and P- 
V characteristics of the P3HT: PCBM polymer solar cell. It is observed that the increase of the temperature 
produces an increase in the short-circuit current (Isc) and a decrease in open circuit voltage (Voc). On the other 
hand, the increase in the temperature leads to a decrease in maximum voltage (Vm) and a slight increase in 
maximum current (Im) and consequently a relative decrease in maximum power (Pmax). Figure 8 (a) shows 
Effect of temperature on the I-V. Figure 8 (b) shows Effect of temperature on the P-V characteristics. 








0 0.1 0.2 0.3 0.4 0.5 0.6 
Voltage (V) Voltage (V) 


(a) (b) 


Figure 8. Effect of temperature on the (a) I-V and (b) P-V characteristics 


Figure 9 represents the variation of the gap energy and the energy conversion efficiency of PSC as a 
function of temperature. It is seen that both the cell efficiency n(%) and the bandgap energy E, decreases 
when the temperature T increases as described by Varshni’s empirical expression [42]. Figure 9 (a) shows 
Efficiency. Figure 9 (b) shows band gap energy of PSC as a function of temperature. 
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Figure 9. (a) Efficiency and (b) band gap energy of PSC as a function of temperature 


aT? 


E (T) = E0 -55 


(17) 
Where E; (0) is the bandgap energy at theoretical zero temperature of material, a and B are material 


constants. As a suggestion, to enhance the PSC’s performance could be achieved through the choice of new 
material in the solar cell domain like perovskite material which is a hybrid of organic and inorganic material 
[43], [44]. 
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6.2. Series resistance effect 

The effect of the series resistance R, on the I-V and P-V curves of the solar cell are illustrated in 
Figure 10. We remark when the series resistance Rs increases, the short-circuit current /,, 1s reduced but the 
open-circuit voltage Vo is not changed. This influence is characterized by a decrease of the slope in the area 
of high voltages. In fact, the series resistance R, characterizes the losses of proper semi-conductor resistance 
and the electrode-organic solar cell contacts ( Al/OSC). As an effect, a high of this resistance decrease 
appreciably the voltage and the output current, and consequently will limit the energy conversion efficiency. 
Besides, when the resistance arises, we perceive a decrease of the output power of the PSC in a considerable 
manner. Figure 11 (a) shows Effect of series resistances on the I-V. Figure 13 (b) shows Effect of series 
resistances on the P-V characteristics 


Current (A) 
Power (W) 





0 0.1 0.2 0.3 0.4 0.5 0.6 
Voltage (V) Voltage (V) 


(a) (b) 


Figure 10. Effect of series resistances on the (a) I-V and (b) P-V characteristics 


6.3. Shunt resistance effect 

These Figure 11 (a) shows illustrates the effect of the shunt resistance Rsp on the I-V. Figure 11 (b) 
shows Effect of shunt on the P-V characteristics of PSHT: PCBM solar cell. We notice that when Rsṣp arises, 
the short-circuit current Ie decreases while the open-circuit voltage V,. remain constant; moreover, this 
influence is characterized by an increase of the slope in the low voltages region. Besides this, the output 
power increases when Rsp decreases. Indeed, the shunt resistance of the D/A junction causes the leakage 
phenomenon and then leads to the losing a part of the photocurrent Ipp in the inner of solar cell device. As a 
result, when the parallel resistance is as high as possible, the maximum power conversion efficiency is 
decreased. 
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Figure 11. Effect of shunt resistances on the (a) I-V and b) P-V characteristics 
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6.4. Influence of the ideality factor 

These Figure 12 (a) shows the influence of the ideality factor n on the I-V and Figure 12 (b) shows 
the influence of the ideality factor n on the P-V characteristics of the solar cell under the Standard Test 
Conditions. It can be seen that the decrease of the ideality factor n causes a reduction in the open circuit 
voltage Voc, while the short-circuit current Isc remains constant. 
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Figure 12. Effect of ideality factor n on the (a) I-V and (b) P-V characteristics 


6.5. Illumination effect 

Figure 13 represents the current-voltage and power-voltage characteristics of a polymer solar cell 
for different solar irradiation ranging from E = 200W/m? to E = 1400 W /m?. We notice that the short- 
circuit current /,- increases significantly comparing to the open circuit voltage Ve. In fact, the short-circuit 
current varies linearly as a function of the illumination, whereas the open circuit voltage varies 
logarithmically [29], [30] as a function of the illumination. As a result, the increase of solar irradiation causes 
an increase in the maximum voltage Vm, the maximum current Im and subsequently a net increase in the 
maximum power Pmax. Figure 13 (a) shows Effect of illumination on the I-V. Figure 13 (b) shows Effect of 
illumination on the P-V characteristics. 
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Figure 13. Effect of illumination on the (a) I-V and (b) P-V characteristics 


7. CONCLUSION 

In this paper, we have studied the dependence both of the intrinsic electrical and the external 
parameters (illumination, temperature, series and parallel resistances, and ideality factor) on the performance 
of polymer solar cells (PSCs) based on the poly(3-hexylthiophene) (P3HT) and [6,6]-phenyl-C61-butyric 
acid methyl ester (PCBM). The single diode model is used to represent the PSC, which is identified by the 
five electrical parameters (Ipn, Io, Rs, Rsh, and n). The proposed method performed by the combination of least 
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mean square method (LMS) and the Newton Raphson’s method to determine the final optimal values of the 
PSC’s parameters. 

The effect of temperature on both the I-V and P-V characteristics, energy band gap, and the energy 
conversion efficiency at different temperatures from 10°C and 90°C has studied, and the influence of the 
electrical and environmental parameters on both of the I-V and P-V characteristic’s PSC. The obtained 
results show that the calculated currents are in a good agreement with the experimental data, and the 
difference between both of the calculated and the experimental currents does not exceed 10% referring to the 
RMSE statistical error. Concerning the effect of the temperature, we have found an increase of Isc, a decrease 
of Voc with increasing of the temperature. However, the variation of the intrinsic electrical parameters shows 
a dependence of the PSC’s efficiency, that we can use in the improvement of PSC performance by taking the 
optimal intrinsic and external parameters. 
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